Localized Structures Embedded in the Eigenfunctions of Chaotic Hamiltonian Systems 
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We study quantum localization phenomena in chaotic sys- 
tems with a parameter. The parametric motion of energy 
levels proceeds without crossing any other and the defined 
avoided crossings quantify the interaction between states. We 
propose the elimination of avoided crossings as the natural 
mechanism to uncover localized structures. We describe an 
efficient method for the elimination of avoided crossings in 
chaotic billiards and apply it to the stadium billiard. We 
find many scars of short periodic orbits revealing the skeleton 
on which quantum mechanics is built. Moreover, we have ob- 
served strong interaction between similar localized structures. 
PACS numbers: 05.45. +b, 03.65. Sq, 03.20.+i 

The connection between individual eigenstates and 
classical invariant sets in chaotic systems is far from be- 
ing fully understood at present. A first step was given by 
Berry |jj and Voros M twenty years ago. They conjec- 
tured that the eigenstates are associated with the whole 
energy surface explored ergodically by the orbits and 
then, the eigenstates would be locally like random super- 
positions of plane waves. That picture is supported by 
Shnirelman's theorem || and the pioneering numerical 
work by McDonald and Kaufman [Q on the Bunimovich 
stadium billiard - an ergodic system |J. However, the 
paper by Heller j^] modified that point of view. He em- 
phasized that a large number of high excited eigenfunc- 
tions for the stadium billiard has an enhancement along 
the shorter periodic trajectories (scars). In 1988, Bogo- 
molny Q] developed the semiclassical theory of scars and 
explained the extra density near unstable periodic trajec- 
tories. He represented the probability density as a sum 
over a finite number of periodic orbits (see also J8|). 

Another unexpected localization phenomenon was ob- 
served in the stadium billiard for the first time; the fa- 
mous bouncing ball (b.b) states ||. They correspond 
to the family of neutral periodic orbits living in the sta- 
dium. Localization is stronger than in the case of scars 
and survives the semiclassical limit p0| . The semiclassi- 
cal theory of b.b was developed recently by Tanner [ |Ll| . 

It is easy to observe b.b states in the stadium billiard, 
but it is practically impossible to find complete series 
of b.b states. On the other hand, it is very difficult to 
find a family of scars. That is, given a scar with wave 
number k which is associated to a particular periodic 
orbit of lenght C, we expect to see a family of scars of 
the same orbit at k + 2-nnJ L with n an integer. But 
in general, a sequence of similar scars is not observed. 
Such a phenomenon appears only for some short periodic 
orbits in some regions of the spectrum but it is impossible 



to define a systematic of scars |12j . 

We believe that in order to study the mechanisms giv- 
ing rise to localization it is necessary first to clean the 
spectrum. The individual properties of two well defined 
states depending on a parameter are mixed when they 
collide in an avoided crossing |jq ]. Then, we expect that 
the elimination of avoided crossings (a.c) in systems gov- 
erned by a parameter dependent Hamiltonian provides 
a mechanism for localization. But in a chaotic system 
it is impossible to remove all a.c because an individual 
eigenstate interacts strongly with other states and the 
interaction is of long range. However, we will extract 
only the a.c with highest curvatures. Recently, it was 
stressed that the distribution of highest curvatures 
has a universal character and the tail of low curvatures 
is a particular property of the system possibly associated 
with scarring phenomena. Precisely, we want to remove 
the universal property of chaotic eigenfunctions related 
to the Shnirelman's theorem, and to retain fluctuations 
revealed like phase space localization. On the same foot- 
ing, Takami jl5) has suggested that a.c are contributed 
by long periodic orbits. In the light of this, by removing 
a.c, scars of short periodic orbits would arise. 

We have developed a simple method to remove a.c in 
chaotic billiards governed by a shape parameter. Let's 
<^t(r) and 0„(r) be normalized eigenfunctions of a pla- 
nar billiard with Dirichlct conditions on the boundary C. 
In ref. Jl6[ ] the following quasiorthogonality relation was 
shown 
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where fc M and k v are the corresponding eigenwave num- 
bers, r n = r.n with n the outward normal unit vec- 
tor to C, and s is arc lenght round C. Relation (|l|) 
defines an effective Hilbert space whose dimension is 
of order k, and where an eigenstate is represented by 
ip(s) = dcj)/dn(r(s)) ; the eigenfunction in this space. 

Now we deform the boundary by changing a shape pa- 
rameter £ (£ = £o for C). Using equation (|l|) it is possible 
to show |l7j that eigenenergies (fc 2 ) and eigenfunctions at 
different values of the parameter are connected through 
the following parameter dependent Hamiltonian written 
in the basis {ffi} of eigenfunctions at £q, 
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K v = -(jf r' n ^ y v ds) x [1 + (k^-K)O(k- 1 )} (3) 

where primes indicate derivation with respect to I. Re- 
lation (||) was obtained previously [Q for the case when 
k a = k lJ . By diagonalizing H{£\) we obtain the eigenen- 
ergies fc^ and the eigenfunctions (p a at Suppose that 
at l\ the states [i and f collide in an a.c. In ref. (n]] an 
efficient way of determining a.c was described in terms of 
the coefficients C^(£) =< (p a \d(p v /d£ > which define a 
driven evolution of the system. Using equation (j|) and 
perturbation theory we obtain C nv {£) = H'/(k? l — kl) , 
where \i and v are associated to eigenstates of H{£) 
(the adiabatic basis). Around £\, behaves like the 
Lorentzian function (£ mt /2)/[£ 2 nt + (£ - £i) 2 ] , with 
tint — l/2C u „(£i) . The area in the energy spec- 
trum, where the a.c is relevant, can be estimated by 
— Ak 2 x £ int . Then, if the area is lower than a 
prescribed value v 

V - (*2 - k 2 v )VK»\ <v = rnl ln2 4 /^ , (4) 

we eliminate the a.c by transforming the diagonal matriz 
of #(4) as follows: k 2 (new) = kl(new) = (k' 2 , + kl)/2 . 
All a.c which satisfy (Q) are eliminated in the increasing 
order of their area. Criterion (Q) has also a simple dy- 
namic interpretation. If the billiard contains a particle 
of mass m and the boundary is driven with a velocity £, 
the a.c is eliminated when the Landau-Zener probability 
transition is greater than one half. 

The new Hamiltonian is less chaotic than the previ- 
ous one, and the adiabatic basis of the new Hamilto- 
nian is more adequate to study nonadiabatic effects p0| . 
For nearly integrable billiards, this new adiabatic basis 
goes to the diabatic one of the original Hamiltonian. For 
chaotic billiards, the elimination of a.c between nearest 
neighboring levels gives rise to a.c between distant levels 
(see Fig. |(a)). 

We have studied the desymmetrized stadium billiard 
with radius r and straight line of length a @. The 
boundary only depends on the shape parameter £ = a/r 
(the area is fixed to the value I+tt/4). Fig. [l] compares 
the approximated spectrum obtained from eq. (||) 
with the exact one; the agreement is excellent. We find 
structures surviving parametric variation (straight lines 
interrupted by small a.c). The most evident being b.b 
states, e.g. states 6, 14 and 18 of Fig. |](a). But in gen- 
eral any localized eigenfunction has the same behavior. 
In Fig. ^(a) we show the spectrum of the transformed 
Hamiltonian by elimination of all a.c satisfying eq.(|^), 
with v = 1.5. At this stage all b.b structures embedded 
in the spectrum emerge clearly, b.b states are identified 
by the quantum numbers {n x ,n y ) counting excitations 
in the horizontal and vertical direction respectively; n y 
labels b.b series. The states 14, 18 and 23 of Fig. |(b) 
correspond to (n x , 15), with n x = 1,2 and 3 resp., and 
states 6, 16 and 28 to (n x , 14), with n, x — 6, 7 and 8. At 



this stage also many scars of short periodic orbits clearly 
appear. They are the states 1, 2, 4, 7, 8, 9, 12, 17, 20, 21, 
22, 26, 28, 29 and 30 of Fig.J|(b); the associated periodic 
orbits are identified in Fig. |[ In Fig. ||(b) we have elim- 
inated avoided crossings up to v — 5 giving rise to new 
localized structures; states 3, 5, 10, 11, 13, 24 and 27 of 
Fig. ||(b). We have not identified the states 15 and 25 
(probably contributed by the whispering gallery family 
of periodic orbits) and the state 19. 

The structures in Fig. |3|(b) are obtained from the spec- 
trum (Fig. |^(a)) by an orthogonal transformation which 
reduces the parametric interaction among states; that is 
, the non-diagonal elements of the deformation matrix 
H' are reduced considerably in the new basis. The new 
states are characterized by the mean energy < k 2 > and 
the dispersion a measured in units of the mean energy 
levels. If (7 is lower than one, the new state has high 
probability of appearing in the spectrum. For a = 1 the 
associated structure lives a time equal to the Heisenberg 
time. For scars of periodic trajectories we expect a to 
increase with k according to the Shnirclman's theorem. 
For b.b states with fixed n x , we expect a to go to zero 
when n y go to infinite. 

The non-diagonal elements of the original Hamiltonian 
in the new basis give us the interaction between localized 
structures. We have observed strong interaction between 
similar structures. For example, b.b states (1,15) and 
(2, 15) have low interaction with the other states (a.c are 
very small -see Fig. |l|), except with the b.b state (7, 14) 
and with the states 10 and 24 of Fig. ||(b). Strong inter- 
action was also observed between scars of short periodic 
orbits which are close in phase space. For example, states 
4 and 22 have strong interaction with states 2 and 17 re- 
spectively, and states 10 and 24 with states 11 and 27 
respectively. In some situations, this strong interaction 
probably decides the occurrence of a defined structure. 
For example, states 4 and 22 are built principally with 
orbit (j), but the Bohr-Sommcrfcld quantization rule for 
this orbit (fc — 0.6999 (n + 1/4)) predicts a scar at k ~ 
47.07 which is not observed. However, the second contri- 
bution, given by orbit (d) (fc = 1.3013 (n + 1/2)), does 
not predict such scar. In fact, in the range 35 < fc < 55 
we have observed structures like 4 and 22 only for those 
values of k which satisfy approximatelly quantization for 
(j) and (d) simultaneously. A deep study of this beat- 
ing phenomenon and the interaction between periodic or- 
bits, which appear very strong and can introduce new 
insight into the semiclassical theory of chaotic systems, 
is presently under way ML 

We transform the eigenfunctions in order to reduce the 
parametric interaction between states. From a semiclas- 
sical point of view, the new structures are not contributed 
by long periodic orbits. The study of the new states in 
terms of short periodic orbits is amenable to the tech- 
nique developed recently in ref. flpJ2| ]. 

Finally, we mention that the large parametric corre- 
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lation observed by Tomsovic |22| in the stadium billiard 
(random matrix theory predicts no correlations) is simply 
understood in terms of these localized structures. On the 
other hand, diffusion and dissipation in mesoscopic sys- 
tems are affected by these localized structures and the 
theory developped by Wilkinson |23j based on random 
matrix theory needs to be revised. 

We would like to thank M. Saraceno for useful sugges- 
tions. E.V. acknowledges the hospitality of U. Smilansky 
and The Wcizmann Institute of Science where part of this 
work has been done. 



FIG. 2. (a) The spectrum of the tranformed Hamiltonian 
by elimination of all avoided crossing satisfying eq. with 
v — 1.5. The scarred states corresponding to Fig.^| (b) are 
shown, (b) Idem (a) for v — 5. 

FIG. 3. (a) Linear density plots of the eigenfunctions for 
the stadium billiar with I — 1 and area l + 7r/4. The numbers 
below each plot are the label (left) and the wave number k. (b) 
Linear density plots of the eigenfunctions for the transformed 
Hamiltonian by elimination of all avoided crossing satisfying 
eq.(Q), with v = 5. Letters on the top identify periodic orbits 
of Fig. |J with h igher contribution to each localized state. The 
label, V< k 2 > and the dispersion a are showm below each 
plot, a is measured in units of the mean energy levels. 
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FIG. 4. Several short periodic orbits of the desymmetrized 
stadium billiard with 1=1. 



FIG. 1. The approximated spectrum (solid lines) obtained 
from eq. (2) is compared with the exact spectrum (dots). 
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